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Quantum phase transitions in the one-dimensional extended quantum compass model in transverse 
field are studied by using the Jordan- Wigner transformation. This model is always gapful except 
at the critical surfaces where the energy gap disappears. We obtain the analytic expressions of all 
critical fields which drive quantum phase transitions. This model shows a rich phase diagram which 
includes spin-flop, strip antiferromagnetic and saturate ferromagnetic phases in addition to the phase 
with anti parallel ordering of spin y component on odd bonds. However we study the universality 
and scaling properties of the transverse susceptibility and nearest-neighbor correlation functions 
derivatives in different regions to confirm the results obtained using the energy gap analysis. 



PACS numbers: 75.10.Jm 

I. INTRODUCTION 

In the last decade, through the extensive experimental 
and theoretical works, the role of orbital degree of free- 
dom in determining the magnetic and transport prop- 
erties of transition-metal oxide materials has been rec- 
ognized extensively^^. The complex interplay among 
the orbital, charge, spin, and lattice degrees of freedom 
makes their phase diagrams extremely rich and leads var- 
ious fascinating physical phenomena. For instance, fer- 
roelectricity, colossal magnetoresistance, and charge or- 
dering are the results of the orbital degeneracy in d-shell 
transition metal oxideiPl. 

A simplified model which described the nature of the 
orbital states in the case of a twofold degeneracy is the 
Quantum Compass Model (QCM) 5 . First, the model 
has been used to describe the Mott insulators with or- 
bit degeneracies. It depends on the lattice geometry 
and belongs to the low energy Hamiltonian originated 
from the magnetic interactions in Mott-Hubbard systems 
with the strong spin-orbit coupling^. In QCM the or- 
bital degrees of freedom are represented by pseudospin 
operators and coupled anisotropically in such a way as 
to mimic the competition between orbital orderings in 
different directions. For simplicity, the one-dimensional 
(ID) QCM, is constructed by antiferromagnetic order of 
X and Y ps eud ospin components on odd and even bonds, 
respectively 7 8 . In addition, the ID QCM is exactly the 
same as the ID reduced Kitaev model--. Brzezicki et al. 
obtained an exact solution for the ground-state energy, 
reveals that the ID QCM exhibits a first-order transition 
between two disordered phases with opposite signs of cer- 
tain local spin correlators. Intriguingly, this first-order 
transition was found to be accompanied by a diverging 
correlation length for spin correlations on one sublattice 7 . 

Moreover, the extended version of the ID QCM, ob- 
tained by introducing one more tunable parameter, has 
been studied by Eriksson et alW^. They have identified 
four distinct ground state phases, are separated by two 



intersecting transition lines. One of them defines a line 
of second-order Ising-like transitions, while the other is a 
line of first-order transitions. The point of intersection, 
where the first-order quantum phase transition identi- 
fied by Brzezicki et al takes place, defines a multicritical 
point. They show that diverging of correlation length 
for certain spin correlations, has recognized a natural ex- 
planation of the multicriticality of the transition point. 
However, their results for the entanglement show that 
the only effect on the ground state when going through 
the first order transitions is that a correlation function 
for neighboring spins on odd bonds changes sign, with- 
out any effect on the entanglement measures. First-order 
quantum phase transitions (QPT) are generally associ- 
ated with a discontinuity in concurrence, but accidental 
exceptions to this rule are possible. To the best of our 
knowledge, the QCM in a transverse field has not been 
studied so far, except on the first order transition lineal 
They have shown that the energy gap does not disap- 
pear in the presence of the transverse-field even in the 
thermodynamic limit. 

In this paper, we study the ID extended quantum com- 
pass model (EQCM) in a transverse magnetic field. The 
exact solution is obtained by using Jordan- Wigner (JW) 
transformation. We show that this model reveal a rich 
phase diagram which includes quantum critical surfaces 
depending on exchange couplings. (We have been made 
aware that similar work is being performed by M. Mo- 
tamedifar, S. Mahdavifar and S. Shayesteh Farjami using 
the exact diagonalization method, personal communica- 
tion.) Moreover, as we have shown in our resent worl^l, 
because of nice scaling properties of correlation functions 
and transverse susceptibility (TS), phase transition can 
be captured from small systems with considerable accu- 
racy without pre-assumed order parameters even for the 
cases where the pairwise entanglement is absent. How- 
ever we have exhibited that the divergence and scaling 
properties of two-body entanglement could be obtained 
by studying the correlation functions properties without 
direct calculation of the entanglement. So we will study 
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the nearest neighbor correlation (NNC) functions and TS 
of this model near a quantum critical points (QCP). 
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II. HAMILTONIAN AND EXACT SOLUTION 

Consider the Hamiltonian 

N' 
n=l 

+ Hvln-l+vL)}- (1) 



a = 2ft, 2 + J 2 + jf + L 2 + 2Lx J 2 cos fc, 

6 = 4[( J x J 2 - h 2 ) 2 + 2 Jii 2 ( ^ J 2 - h 2 ) cos fc + J 2 L 2 ]. 

The ground state (-Eg) an< i the lowest excited state 
(Ee) energies are obtained from Eq.(j3|, 

E G = -W(El + EH E E = -\Y,(El-Ei), 



where J\ and J 2 are the odd bonds exchange cou- 
plings, L\ is the even bond exchange coupling and 
N = 2N' is the number of spins. We assume periodic 
boundary conditions. The above Hamiltonian (Eq. ([IJ) 
can be exactly diagonalized by standard Jordan- Wigner 
transformation^ as defined below, 



b t 



b t 



2b t b J 



b J = * 



which transforms spins into fermion operators Cj . 

The crucial step is to define independent Majorana 
fermions^ at site n, = c 2n _i and c£ = c 2rl . This 
can be regarded as quasiparticles' spin or as splitting the 
chain into bi-atomic elementary cells^. 

Substituting for crj, <j| and erj (j = 2n, 2n— 1) in terms 
of Majorana fermions with antiperiodic boundary condi- 
tion (subspace with even number of fermions) followed 
by a Fourier transformation, Hamiltonian Eq. ([I]) (apart 
from additive constant), can be written as 



H+ = 



2h(cfcl + $4) + h.c. ,(2) 



where J = 
and k = ±jp 



{Jx - J 2 ) - L ie ik , 
, 0' = 1,3, ••• ,N'- 



L = 
-!)• 



Ji + Ja) + L ie l 



It should be pointed out that although the GS in pe- 
riodic and antiperiodic boundary conditions are slightly 
different in the finite-size system, they are identical in the 
thermodynamic limit and the essential features in finite 
size are also not altered qualitatively. 

Finally, diagonalization is completed by a four- 
dimensional Bogoliubov transformation connecting 



quasiparticles, 
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c u and obtain two different kind of 



H = 



E 

k 



Eliik 



at 



(3) 



where could be written as a function of a and 6, 



Er 



= -2^V a + A E E = -2j2\a-Vb (4) 



k>0 



k>0 



It is clear the ground state is separated from the lowest 
energy pseudospin excitation by a pseudospin gap A = 
\Ee — Eq |, which vanishes at h = W Ji( J 2 + L{) and 
fin = v/ J\{ J 2 — L\) in the thermodynamic limit. 

It should be stressed here that the exact spectrum and 
the pseudospin gap are the same as that obtained by 
Eriksson et al. for h = using a different method^. 

So, the quantum phase transition (QPT) which could 
be driven by the transverse-field, depending on exchange 
couplings, occurs at /iq and h w . 



III. PHASE DIAGRAM 

The complete phase diagram of the extended compass 
model without transverse magnetic field has been re- 
ported in Refs. [TU] and [TS]. They have shown that the 
first-order transition occurs at multicritical point where 
a line of first-order transition (Ji/L\ = 0) meets with a 
line of second order transition (J 2 /Li = 1). Also, There 
are four gapped phases in the exchange couplings' space, 

• (I) Ji > 0, J 2 < 1: In this region the ground state 
is in the Neel phase along the x axis. 

• (II) Ji > 0, J 2 > 1: In this case there is antiparallel 
ordering of spin y component on odd bonds. 

• (III) Ji < 0, J 2 > 1: In this case there is parallel 
ordering of spin y component on odd bonds. 

• (IV) Ji < 0, J 2 < 1: In this region the ground 
state is in the strip antiferromagnetic (SAF) phase. 

Phase diagram of extended quantum compass model in 
transverse field has been shown in Fig. ([I]). Depending 
on exchange couplings, the transverse field could drives 
the phase transition at ho and h„ where the energy gap 
vanishes (For simplicity we take Li = 1). 
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FIG. 1: (Color online) Phase diagram of the extended com- 
pass model in the transverse magnetic field. For Ji > 0, J2 < 
1, the front and top side of the purple convex surface is spin- 
flop phase (I) and the back side ( Ji > 0, J2 > 1) specified 
by antiparallel ordered of spin y component (II). The blue 
checkerboard pattern represents the boundary between spin- 
flop phase (I) and saturate ferromagnetic phase (III). In the 
case of Ji < 0, J2 < 1 there are two phases, the strip antifcr- 
romagnetic phase (IV) which exists below the green convex 
surface and the saturate ferromagnetic phase (III) which is 
above it. For Ji < 0, J2 > 1 saturate ferromagnetic phase 
(III) attends the phase diagram. The first-order transition 
appears just at Ji = 0, h — line. 



Fig. (pi) shows the absolute value of transverse magne- 
tization (TM) and NNC functions on odd and even bonds 
for infinite system size in region (I) ( J\ = 1, J2 = 0.8). In 
this region tuning the magnetic field dictated the system 
to fall into saturated ferromagnetic (SF) phase. The spin 
flop-SF phase transition occurs at h c = ho (blue checker- 
board curved plane in Fig. (Ji}) which under this surface 
ground state is in the spin flop phase (the Need ordered 
along the axis where is perpendicular to magnetic field 
is called spin flop). It is seen in Fig. ^ that the on- 
set of magnetic field sets up the TM (M z ) immediately 
and continuously increases with increase in h to satu- 
rate value (|M Z | = 1). However, the antiparallel ordered 
of spin x and y components on odd (G* x , G^ v ) and even 
(Gff ) bonds reduces by increasing the magnetic field and 
goes to zero for h — > 00. 

In region (II) the gap decreases by increasing the mag- 
netic field and goes to zero at the lower critical field 
h ci = h„ (purple checkerboard curved plane in Fig. ([!)) 
and beyond this critical field the energy gap immediately 
appears with increase of the magnetic field. This process 
continues until the upper critical field h C2 — ho (blue 
checkerboard curved plane in Fig. (JlJ) at which the en- 
ergy gap vanishes and becomes once again gapped upon 




FIG. 2: (Color online) The transverse magnetization and dif- 
ferent components of nearest-neighbor spin correlation func- 
tions on even and odd bond for Ji = 1, J 2 = 0.8. 



the h C2 . Fig. (|3| shows the TM and NNC functions ver- 
sus the magnetic field in the region (II) (J x = 1, J 2 = 2). 
It manifests that under the lower critical field (h < h Cl ) 
the antiparallel ordered of spin x and y components on 
odd and even bonds stay quite unchanged. However the 
TM is zero for h < h Cl . So the ground state's antiparallel 
ordering of spin y component remains unchanged under 
h ci . Beyond the h Cl the TM and NNC functions undergo 
a strong qualitative change and antiparallel ordered on 
odd and even bonds tends to zero as the magnetic field 
increases. Increasing the magnetic field saturates the TM 
and disappears the antiparallel ordered spin x component 
on odd bound at h C2 = ho , while the antiparallel ordered 
spin y component on odd bound and antiparallel ordered 
spin x component on even bond have a nonzero values 




FIG. 3: (Color online) The TM and NNC functions of spin 
components for Ji = 1, J2 = 2. 
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FIG. 4: (Color online) The TM and different components of 
nearest-neighbor spin correlation functions on even and odd 
bond in region (III) (Ji = —1, J2 = 2). 



and tend to zero for h — > 00. 

A surprising result occurs in the intermediate region of 
the magnetic fields h Cl < h < h C2 where increasing the 
magnetic field enhances the antiparallel ordered spin x 
component on even bond up to a maximum and then de- 
creases gradually, while decreases the other antiparallel 
ordered. So we predict that the gapped spin-flop phase 
exists in the intermediate values of the transverse mag- 
netic field h Cl < h < h C2 . In other words, in the region 
(II), the magnetic field destroys the ground state's an- 
tiparallel ordering of spin y component on even bond at 
h Cl and sticks the system in the spin-flop phase upon the 
h ci . The spin-flop-SF transition occurs beyond h C2 . 

The TM and NNC functions have been depicted in Fig. 
@ for region (III) {J\ = — 1, J2 = 2). In this case the 
system is in the SF phase and phase transition dose not 
occur by tuning the magnetic field even in the strong 
magnetic field. 

Fig. |pl) shows the TM and NNC functions in the re- 
gion (IV) (Ji = — 1, J-i = 0.8). This region include two 
gapped phases, SAF and SF where separated from one 
another at the critical point h c = h K (green checkerboard 
curved plane in Fig. Q). 

The three-dimensional panorama of TM with respect 
to J 2 and h has been plotted in Fig. ^ for J x = 1. The 
two critical line (ho(J2,h) and h^{J2,h)) at which the 
energy gap vanishes can be described by two assumed 
lines on the two convex part of the surface in Fig. ([5]). 



IV. UNIVERSALITY AND SCALING OF 
CORRELATION FUNCTIONS 

The nonanalytic behavior in some physical quantity is 
a feature of second-order quantum phase transition. It is 
also accompanied by a scaling behavior since the correla- 




FIG. 5: (Color online) The TM and NNC functions of spin 
components for J\ = —1, J2 = 0.8. 
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FIG. 6: (Color online) The three-dimensional panorama of 
TM for Ji = 1. 



tion length diverges and there is no characteristic length 
scale in the system at the critical point. As we previously 
mentioned, correlation functions and TS show the uni- 
versality and scaling around the QCP and could capture 
QCP. However, studying the NNC functions behaviors 
could reveal the scaling and universality of entanglement 
near the QCP. So, in this section we will study the be- 
havior of NNC functions derivative with respect to the 
magnetic field and TS to confirm the previous results. 

In Fig. Q the derivative of G v a v with respect to the 
magnetic field has been shown for different system sizes 
in region (I) (Ji = I, J2 = 2). For infinite lattice size 
dG v v /dh diverges as the critical point is touched, while 
there is no divergence for finite lattice sizes. As the size 
of system becomes large, the derivative of NNC functions 
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FIG. 7: (Color online) Evolution of the versus h for 

different system sizes in region (I) for Ji = 1, J 2 = 0.8. 



FIG. 9: (Color online) Finite-size scaling of dG yy /dh for dif- 
ferent lattice sizes. The curves which correspond to different 
system sizes clearly collapse on a single curve. 
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FIG. 8: (Color online) Scaling of the maximum of 
systems of various sizes. Inset: Scaling of the position (h,Max) 
of dG for different-length chains. 



tends to diverge close to the critical point. More informa- 
tion can be obtained when the maximum values of each 
plot and their positions are analyzed. The position of the 
maximum {huax) of dG v v /dh tends toward the critical 
point like h Max = h c - N~ 9 (6 = 1.72 ± 0.03) which has 
been plotted in the inset of Fig. Q. 

Moreover, we have derived the scaling behavior of 
\dG^, v / dh\h Max versus N. This has been plotted in Fig. 
^ which shows a linear behavior of \dG yy /dh\h Max ver- 
sus ln(iV). The scaling behavior is \dG yy / dh\h Max = 
rln(AT) with r = 0.15 ± 0.01. To study the scaling be- 
havior of G yy around the critical point, we perform finite- 
scaling analysis, since the maximum value of derivative of 
G v a v scales logarithmic. According to the scaling ansatz, 
the rescaled derivative of G yy around its maximum value 
hMax is just a function of rescaled driving parameter such 



dG yy dG yy , 



dh 



dh 



F(N^ v {h-h Max )), (5) 



G|f and 



where F(x) is a universal function. The manifestation of 
the finite-size scaling is shown in Fig. (Joj) . It is clear that 
the different curves which are resemblance of various sys- 
tem sizes collapse to a single universal curve. Our result 
shows that v = 1 ± 0.001 is exactly correspond to the 
correlation length exponent of Ising model in transverse 
field [y = 1). 

A similar analysis can be carried on G% 
TS (x z ). Our calculations show that the non-analytic 
and scaling behavior of NN correlation functions are the 
same as TS does. It is important to mention that NNC 
functions and TS show the logarithmic divergence near 
the QCP. Our results is different from the reported result 
in Ref . [TT] . They have found that TS shows a power-law 
behavior close to the QCP with the exponent 7 = 1.78 ± 
0.05. Since h is analogous to temperature in classical 
systems, we expect h\ z to be equivalent to the specific 
heat in the 2D classical Ising model (J2 = 0). So the 
reported exponent in Ref. [TT] belongs to susceptibility 
in the x-direction not to TS. 

We have plotted the derivative of G yy in region (II) 



(Ji = 1, J2 = 2) versus h in Fig. (10 1 for different lat 



tice sizes which shows the singular behavior as the size 
of the system becomes large. As it manifests the diver- 
gences of dG yy /dh occur at h Ct — 1 and h C2 — v3 where 
exactly correspond to the critical points that obtained us- 
ing the energy gap analysis (h Cl — h w , h C2 — ho). A more 

dG yy 

detailed analysis manifest the linear behavior of .? at 
the first maximum point {huaxt) versus ln(iV) where has 
been plotted in Fig. (11). The exponent for this behavior 
is 7i =0.14 ±0.02. 

Moreover, we have shown that the position of the first 
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FIG. 10: (Color online) The first order derivative of dG yv /dh 
as a function of h for various system size in region (II) (Ji = 
l,Ja = 2). 
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FIG. 11: (Color online) The scaling behavior of the first max- 
imum point of dG dh for different-length chain in region (II). 

Inset: Scaling of the position (ftifaii) of for different- 

length chains (Ji = 1, J2 = 2). 



maximum (fiMaxi) of dG v v /dh goes to the first critical 
point, such as hj^j axi — h Cl — N~ 61 with 6*i = 1.59 ±0.01 
(Fig. (11 inset)). The similar investigation shows the 

. dG yy 

scaling behavior of ,? at the second maximum point 
{%Maxi) where has been presented in Fig. \12\ . It spec- 
ifies a linear behavior of dG ,? \% Max versus ln(iV) with 



the same exponent as d ° treat at the first maximum 
point (T2 = 0.14 ±0.02). However, the second maximum 
position of correlation functions and TS show the scaling 
behaviors with the same exponents (Fig. (12), inset) as 
the position of the first maximum of correlation functions 
and TS do. 



FIG. 12: (Color online) The logarithm of the second maxi- 
mum of dG dh I versus the logarithm of chain size, ln(./V), which 
is linear and shows a scaling behavior (Ji = 1, J2 =2). Inset: 
The scaling behavior of hM aX2 in terms of system size (N) 
where hyiaxi is the position of second maximum in Fig. ( 10 1 
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FIG. 13: (Color online) The finite-size scaling analysis for 
the case of logarithmic divergence around the first maximum 
point (liMaii) for J\ = 1, J% = 2. The NNC function, consid- 
ered as a function of system size and coupling, collapses on a 
single curve for different lattice sizes. 



around its first and second maximum points in figs. ( 14) 
and ( 15 1 respectively. They show that the NNC functions 



We illustrate the finite size scaling behaviors of 



dG y " 
dh 



can be approximately collapsed to a single curve. These 
results show that all the key ingredients of the finite size 
scaling are present in these cases too. In this cases scaling 
is fulfilled with the critical exponent v — 1 in agreement 
with the previous results and universality hypothesis. A 
similar analysis show that the universality and scaling 
behavior of G% x , G|f and TS (x z ) are the same as each 
other in region (II). 

In Fig. ( 15 ) the derivative of G X X with respect to mag- 
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FIG. 14: (Color online) A manifestation of finite-size scal- 
ing of around the second maximum point (ftjnoi 2 ) f° r 
various system sizes in region (II). 
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FIG. 15: (Color online) The derivative of dG vy versus h in 
region (III) for J\ = — 1, J% = 2. Even in the thermodynamic 
limit no singularity is observed. 



netic field has been shown versus h for J\ = — 1, J% = 2 



(region (III)). As it is clear the correlation functions does 
not show any singularity in this region even at the ther- 
modynamic limit. This is justify our previous finding 
that shows no transition in this region. We have also 
investigate the TS and NNC functions behavior in the 
region (II) . Our examination show the logarithmic diver- 
gence and scaling behavior of them close to the critical 
point where has been obtained using the energy gap anal- 
ysis. 

V. SUMMARY AND CONCLUSIONS 

In this work we have studied the quantum phase tran- 
sition in the one-dimensional extended quantum compass 
model in the presence of a transverse magnetic field. We 
have shown that there are strip-antiferromagnetic, spin- 
flop and saturate ferromagnetic phases in addition to the 
phase with anti parallel ordering of spin y component on 
odd bonds where phases separated from each other by 
the critical surfaces. We obtain the analytic expressions 
of all critical fields for the field-induced quantum phase 
transitions (QPT). However we have investigate the uni- 
versality and scaling properties of the nearest neighbor 
correlation functions and transverse susceptibility to con- 
firm the results were obtained using the energy gap anal- 
ysis. The results show that the transverse susceptibility 
and derivatives of the nearest neighbor correlation func- 
tions diverge close to the critical point and exhibit beau- 
tiful scaling law. So we predict that as the correlation 
length diverges at the critical point for an infinite system 
size, the derivative of correlation functions between far 
neighbors could capture the quantum critical point too. 
The obtained exponents (v = 1) and universality behav- 
iors (logarithmic divergence) are nearly the same as those 
in the ID transverse-field Ising modeP^, suggesting that 
these two models share the same universality class. 

Further investigations including blocks and multi-body 
correlations functions and the effect of temperature may 
be interesting to establish a precise comparison between 
universality and scaling behaviors of correlation functions 
and entanglement at critical points. Moreover, dynamics 
of correlation functions may also provide an interesting 
scenario for the discussion of the properties of the correla- 
tion functions and its implications for phase transitions. 
Such topics are left for a future research. 
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